Homoclinic Structure Controls Chaotic Tunnelling 
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Tunnelling from a chaotic potential well is explained in terms of a set of complex periodic orbits 
which contain information about the real dynamics inside the well as well as the complex dynamics 
under the confining barrier. These orbits are associated with trajectories which are homoclinic to a 
real trajectory emerging from the optimal tunnelling path. The theory is verified by considering a 
model double-well problem. 
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Chaos is ubiquitous in nature. For this reason, its effect 
on quantum systems is very well studied [Q and powerful 
semiclassical techniques have been developed to treat it. 
However, tunnelling has largely remained impervious to 
such analysis. We show that a simple principle underlies 
the detailed behaviour of tunnelling in chaotic systems 
- state to state fluctuations in the tunnelling rate are 
determined by a restricted set of complex classical trajec- 
tories which can be constructed straightforwardly using 
known techniques of dynamical systems. 

Previous analyses of chaotic tunnelling exist in con- 
texts where a quantum state is initially localised in an in- 
tegrable region of phase space and tunnels into or through 
a chaotic region Our goal is different; we con- 

sider systems that are almost completely chaotic and 
tunnelling is through an energetic barrier. We then con- 
struct a periodic-orbit theory of tunnelling analogous to 
Gutzwiller's celebrated formula for the density of states. 
It relates tunnelling rates to a set of complex classical tra- 
jectories which can be understood intuitively as follows. 
In any potential barrier, there is an optimal tunnelling 
route determined by a complex trajectory with minimal 
imaginary action. Interference effects in the tunnelling 
rate are produced by trajectories which originate from 
and return to a small region of phase space surrounding 
the minimal trajectory. In chaotic systems, such recur- 
rent trajectories are understood on the basis of "homo- 
clinic intersections". Technical details come later, but 
the main point is that such trajectories are standard in 
dynamical systems theory and are relatively easy to find. 
These orbits dominate tunnelling rates in a manner sim- 
ilar to the way in which they control wavepacket recur- 
rences as discovered in [Q . See also H . 

We treat double well potentials for which V{— x, y) = 
V(x, y) and examine the energy level splittings of quasi- 
degenerate doublets labelled by n (although our general 
conclusions apply equally to resonances in metastable 
wells )■ Detailed numerical calculation will be for 
V(x, y) = {x 2 - l) 4 + x 2 y 2 + [iy 2 with /i = 1/10. It will 
prove useful, instead of dealing directly with the energy- 
level spectrum, to regard q — 1/h as a parameter and ask 



for the values q n for which a given energy E is an eigen- 
value — this leads to calculations with a fixed classical 
dynamics. We then obtain doublets at values q n with 
splittings Aq n between them. We investigate tunnelling 
by constructing the function |l, 



/(<?) = ^2^QnS{q - q n )- 
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The spectrum and tunnelling splittings are encoded in 
f(q). Semiclassically, it can be approximated as a sum 
over pseudoperiodic orbits which start on one side of the 
energy barrier and, after evolution in complex time, fin- 
ish at the symmetric image of the initial point. In Q 
it was shown that when the optimal tunnelling route is 
along a symmetry axis, a quasi-periodic oscillation arises 
in the tunnelling rate which is related to the action of a 
real periodic orbit, also on the symmetry axis. However, 
there remain significant fluctuations which we now ex- 
plain using non- axial orbits homoclinic to this real orbit. 

We can visualise a complex pseudoperiodic orbit as a 
one-dimensional path in complex phase space. However, 
it is not unique — it corresponds to a given contour de- 
scribing its complex time evolution; any deformation of 
this contour gives an equivalent orbit. Because of this 
freedom, a systematic and unambiguous description is 
difficult if we work in full phase space. The situation sim- 
plifies considerably if we restrict dynamics to a surface 
of section. Pseudoperiodic orbits then reduce to single, 
isolated points of intersection with the surface of section 
and classification of them is simpler. An added advan- 
tage is that metastable sytems can be treated within the 
same formalism §,0. 

Let £ denote a surface of section through one of the 
wells, defined in the usual manner of real dynamics (eg. 
x = xo, x > and H = E.) Let F : £ — > £ denote the 
first return surface of section mapping, defined by letting 
trajectories start on £ and evolve under real dynamics 
until they intersect £ once again. Real periodic orbits 
correspond to fixed points of some iterate F r of F. We 
incorporate tunnelling by introducing a second complex 
map T acting on the same section £ as follows. 

Any potential barrier has an orbit 70 which crosses 



it with minimum imaginary action iKq in an imaginary 
time —itq. It evolves in complex phase space except at 
two turning points where it is fully real. There is a real 
point zq on E which evolves to one of these endpoints un- 
der real time evolution. It has a symmetric image <tzq on 
erE which evolves to the other endpoint. (Here a refers to 
the symmetry operation — reflection in x in our exam- 
ple.) The point zq evolves into azo if we integrate along 
a complex time contour C as follows: first the contour 
evolves along the real axis until the trajectory reaches 
the turning point of 70, then it descends parallel to the 
imaginary axis until the turning point on the other side 
of the barrier is reached and finally retraces the evolution 
parallel to the real axis (in negative real time) until azQ 
is reached. The net time of evolution is — itq. 

Wc define the map T by deforming 70. If a point z' 
near z$ is chosen on S, a deformation of C can be found 
such that the final point z" is on cE. By invoking the 
symmetry operation to map erE back to E, we define a 
complex symplectic map T : z' 1— > z = a~ 1 z" from E 
onto itself. The central orbit 70 corresponds to a real 
fixed point z of T while other initial conditions, even if 
real, are mapped to complex images under T . Just as 
the real periodic orbits give real fixed points of F r , the 
complex pseudoperiodic orbits give complex fixed points 
of TF r \ such fixed points define orbits with r oscillations 
in one well before tunnelling across the barrier. 

Following 11 , the classical maps F and T have quan- 
tisations H as unitary and hermitean operators, respec- 
tively, acting on a Hilbert space which quantises E. The 
expression for f(q) as a sum over the fixed points of TF r 
has an interpretation as a sum over traces of such oper- 
ators, as will be shown in a future publication fl(i|| . Ne- 
glecting multiple tunnelling traversals and the uniform 
analysis needed at the bottom and top of the wells, 
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where P r denotes the set of fixed points of TF r . The 
fo(q) term is determined by 70, is monotonic and gives 
the average behaviour of the splittings || . The contribu- 
tions with r > 1 have complex actions S-y with nonzero 
real parts and describe fluctuations superimposed on 
fo(q)- The amplitudes are A 7 = 1 / ^/— det ( Af 7 — 7) 
where M 7 is the complex 2x2 symplectic matrix lin- 
earising TF r about 7 (note that a different phase con- 
vention was used in Q.) The choice of root in A-y can be 
unambiguously assigned using the real dynamics [fl0| . 

In problems with an additional symmetry axis (for us 
it is along y = 0), the point zq is a fixed point of both 
T and F and thus of each TF r . The treatment in ||] 
consisted of treating only this fixed point in each of the 
terms in (^). We now systematically include the other 
fixed points, which are generally complex. 
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FIG. 1. Dashed curve: Fourier transform of g(q) obtained 
quantum mechanically. Upper solid curve: semiclassical pre- 
diction using just the axial orbits. Lower solid curve: theory, 
using the six homoclinic families shown in Fig. 2. Inset: the 
quantum-mechanical Keg(s) and the semiclassical prediction 
in a limited range (note there are two superimposed curves). 
Peaks not accounted for by the lower curve correspond to 
uncomputed secondary intersections. 

We now scale out the smooth exponential variation 
with q by defining g{q) = f(q)/fo(q) which is a function 
of order unity. It is the same as Eq. (|j) but with Aq n 
replaced by x n = Aq n /f (q n ). We find 

00 

g(q) = 1 + ReJ2 E ( 3 ) 

r = l jEPr 

where A'^ is the rescaled amplitude. Note that q only 
appears in the exponents and that the residual action 
S~, — iK(j is predominantly real if Im ~ Kq so that 
g{q) is an oscillatory function. 

In Fig. [I] we show the Fourier transform g(s) of g{q) for 
energy E — 0.9 and q £ [30, 100], chosen to put us well 
within the semiclassical regime. For each r, the fixed 
point of J-F r at zq has action 5 7 = rS*o + iKg, where 
Sq is the action of the on axis real periodic orbit in one 
of the wells. Therefore g(s) exhibits a series of simple 
peaks at the harmonics s = rSo, well reproduced by the 
semiclassical theory. More interesting is the additional 
structure, primarily in the form of a sequence of regularly 
spaced peaks shifted away from the simple harmonics, 
which we now explain using nonaxial orbits. 

Among barrier-crossing orbits, the minimum imagi- 
nary action is that of 70; the contributions of other or- 
bits are suppressed by order exp(— AK/K) where AK is 
the additional imaginary action. To be numerically sig- 
nificant, orbits in the sum Eq. (|3|) should have a small 
AK and so should begin and end near Zq. On the other 
hand, to contain information about the entire spectrum 
as implied in (|l]), they should fully explore the phase 
space. These two requirements are satisfied by trajecto- 
ries homoclinic to zq. Asymptotically, these trajectories 



approach zq under both forward and backward time evo- 
lution as d(t) ~ e _I/ '*', where v is the stability exponent 
of the central orbit and t counts intersections with the 
Poincare section. In practice they are easily calculated 
as the intersections of the stable and unstable manifolds 
of zq (the sets of points in phase space which approach 
Zq under forward and backward time evolution, respec- 
tively.) For a general discussion see IO ]. 

Let (• • • X-2, sc-i> xq, x\, x% ■ • •) be the intersection 
of one such trajectory with S. Then a truncation 
(x-m ■ ■ • X-i, xq, X\, ■ ■ ■ , Xn) is a finite length trajectory 
which is almost a periodic orbit if M and N are large 
because X-m and xn are exponentially close to Zq and 
to each other. By slightly perturbing the initial condi- 
tion, one can find a nearby fixed point of TF r ', where 
r = M + N, by complex Newton iteration. The tun- 
nelling segment of this trajectory is very close to 70 and 
as a result the imaginary part of the action is close to 
the optimal value Kq. In the trace formula, there is then 
a competition of exponentials — the imaginary part of 
the action will decay exponentially with length according 
to Im(S* 7 - iK ) ~ e -»™*{M,N) and this can be sma n er 
than h so that e -i(S-,-iK ) j g n0 nnegligeable. For a fixed 
r, each choice of M (and N — r — M) gives a distinct 
fixed point although only a subset will be numerically sig- 
nificant (the effective number depends on q and is given 
by the criterion that Im(5 7 - iK ) ~ e -™&(M,N) < h y 
While the arguments above are asymptotic in M and N, 
in practice complex orbits could usually be found as long 
as M, N > 2 for the system we examine. 

The near equal spacing of the peaks is because one ex- 
tra iteration of the real map amounts to including one ex- 
tra -F-bounce near zq] to a good approximation this adds 
So to the action. Also note that the homoclinic peaks are 
larger than the peaks at s = rSo beyond about r = 5. 
Each homoclinic peak actually contains the contribution 
of many fixed points, all with approximately the same ac- 
tion. This quasidegeneracy is because there is more than 
one homoclinic trajectory and because of the choice of M 
for a given r as explained above. The amplitude of any 
single orbit initially grows algebraically with length due 
to the increasing quasidegeneracy of the orbit. This is ul- 
timately overcome by the exponential decay with length 
due to the instability of long orbits. In this case the 
maximum is around r = 7. There is additional struc- 
ture which dominates at larger r and which we ascribe to 
secondary homoclinic intersections. This scenario is like 
the calculation of wavefunction recurrences in 0, where 
there is also a selection of trajectories which return to a 
classically small region of phase space. There the region 
is determined by a Gaussian wavepacket, whereas here it 
is by the near-Gaussian kernel of an operator quantising 
T Iqjiotl . Fourier analysis reveals peak structure similarly 
organised around homoclinic orbits . 
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lowed region. W u and W a are the stable and unstable man- 
ifolds whose intersections define the homoclinic trajectories: 
the six distinct ones being labelled a to /. For each one, we 
show the trajectory in configuration space of a correspond- 
ing tunnelling orbit. Actually shown is the full periodic orbit 
which is the double iteration of the pseudoperiodic orbit used 
in calculation. Also, only the real parts are shown; the small 
imaginary components being too small to see. 

We now apply this formalism to the potential V(x, y). 
Shown in Fig. || are the stable and unstable manifolds of 
zq in the x = 0.6 surface of section. Of particular inter- 
est is a sequence of 6 primary intersections labelled a to 
/. (All other intersections shown are either iterates of 
these six or are related by reflection in y and p y .) This 
is a rather rich structure; often there are only 2 primary 
intersections |Tj] . Each of these intersections defines a 
homoclinic trajectory, and a corresponding series of fixed 
points of J r F r . We show the projections in real config- 
uration space of complex trajectories defined for each of 
these intersections by the truncation (M, N) — (3, 5) cor- 
responding to the peak at s w 8.25*0 . 

The symmetries of time-reversal and reflection in y 
lead to degeneracy among the orbits. The families b and 
c are respectively mapped into those of / and e under a 
combination of the two symmetries. Congruent but dis- 
tinct families are obtained if we apply either symmetry 
alone. The family a maps to itself under time-reversal but 
a distinct family is obtained under reflection in y, while 
the converse applies to d. We therefore only have four 
inequivalent families: a, b, c and d with degeneracies 2, 
4, 4 and 2 respectively. This high degeneracy (including 
also the choice M), compared to the single axial orbit, is 
responsible for the relative dominance of the associated 
peaks in Fig. 0. In the resolution presently available, the 
individual contributions of {a, 6, c, d} are not resolved, 



but presumably would be with a wider g-window. 

To get detailed agreement in Fig. [I] we had to include 
a family of "ghost" orbits in addition to the orbits de- 
scribed above. Near the intersections a and b, there is a 
switchback in the stable manifold emanating from / that 
almost produces two additional intersections; at lower 
energies these intersections actually occur and lead to 
two new families of orbits. These can be tracked up to 
E = 0.9, although they become complex. (These orbits 
are analogous to the ghost orbits of (ll|.) One of them 
is exponentially large and is excluded through the Stokes 
phenomenon. The other does contribute and is an impor- 
tant component in the leading edges of the peaks. The 
Fourier transform is then almost completely reproduced 
up to about s = 8Sq. Thereafter additional side peaks 
emerge which can be explained by the secondary inter- 
sections (which would be seen in Fig. |^ if the invariant 
manifolds were extended.) Therefore, a systematic cal- 
culation of homoclinic structure suffices for a complete 
understanding of chaotic tunnelling. 

While the details presented here are limited to systems 
with additional symmetry in y, we claim that the guid- 
ing principle applies generally. If the real extension of 
7o is not periodic, the tunnelling is given by trajectories 
homoclinic to it, though they are not as simply organ- 
ised as here. The important physical notion is that ho- 
moclinic trajectories are asymptotic to the optimal tun- 
nelling route while fully exploring classical phase space. 

There has been significant recent effort to understand 
dynamical tunnelling in which the barriers are not en- 
ergetic but come from dynamical effects such as KAM 
surfaces [^]]3). In particular, an exhaustive semiclassical 
analysis was carried out for one particular system in || . 
The authors considered the propagation of a quantum 
state initially localised in a regular region of phase space 
and used complex trajectories to describe its penetration 
into a chaotic region. Explicit semiclassical formulas were 
developed but it is difficult to make a direct comparison 
to our results. Here we have considered a purely spectral 
quantity in terms of classical invariants — the complex 
periodic orbits — while they considered the propagation 
of a specific quantum state. In terms of the quantities 
considered, namely splittings, the work of [Q is closer 
in spirit. However, because of the complicated, mixed 
nature of the phase space, there are no explicit semiclas- 
sical results of the type derived here. In particular, the 
splittings typically vary over many orders of magnitude, 
even after normalisation by the local mean, unlike what 
we observe. While it would be desireable to develop an 
explicit semiclassical formalism for such systems, it is far 
from obvious how this could be done. 

The problem of the statistics of the normalised split- 
tings will be addressed in a later publication jl(|. For 
now we remark that they are not governed by the Porter- 
Thomas distribution but depend on the specific proper- 
ties of 7o- This is in spite of the fact that splittings are 



formally similar to resonance widths for which one ex- 



pects Porter-Thomas 14 1. This effect is quite general 
and should also apply to certain resonance problems. 



' Present address: Division of Theoretical Mechanics, 
School of Mathematical Sciences, University of Notting- 
ham, NG7 2RD, UK. 
* Unite de recherche des Universites de Paris XI et Paris 
VI associee au CNRS. 

[1] M. C. Gutzwiller, Chaos in Classical and Quantum Me- 
chanics (Springer Verlag, New York, 1990). 

[2] W. A. Lin and L. E. Ballentine, Phys. Rev. Lett. 65, 
2927 (1990); O. Bohigas, S. Tomsovic and D. Ullmo, 
Phys. Rev. Lett. 64, 1479 (1990); Phys. Rep. 223, 45 
(1993); O. Bohigas et al., Nucl. Phys. A560, 197 (1993); 
S. Tomsovic and D. Ullmo, Phys. Rev. E 50, 145 (1994). 
E. Doron and S. D. Frischat, Phys. Rev. Lett. 75, 3661 
(1995); Phy. Rev. E 57, 1421 (1998). 

[3] A. Shudo and K. Ikeda, Prog. Theor. Phys. Supplement 
No. 116, 283 (1994); Phys. Rev. Lett. 74, 682 (1995); 
ibid. 76, 4151 (1996); Physica D 115, 234 (1998). 

[4] S. Tomsovic and E. J. Heller, Phys. Rev. E 47, 282 
(1993). 

[5] A. Ozorio de Almeida, Nonlinearity 2, 519 (1989). 
[6] S. C. Creagh and N. D. Whelan, Ann. Phys. 272, 196 
(1999). 

[7] M. W. Beims, V. Kondratovich and J. B. Delos, Phys. 

Rev. Lett. 81, 4537 (1998). 
[8] S. C. Creagh and N. D. Whelan, Phys. Rev. Lett. 77, 

4975 (1996). 
[9] E. B. Bogomolny, Nonlinearity 5, 805 (1992). 
[10] S. C. Creagh and N. D. Whelan, to appear. 
[11] S. Wiggins, Chaotic Transport in Dynamical Systems, 

Springer- Verlag, New York (1992). 
[12] K. Hirai and E. J. Heller, Phys. Rev. Lett. 79, 1249 
(1997). 

[13] M. Kus, F. Haake and D. Delande, Phys. Rev. Lett. 71 
2167 (1993). 

[14] C. E. Porter (ed.), Statistical Theory of Spectra, Aca- 
demic Press, New York (1965). 



